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Intrinsic, two-dimensional ferromagnetic semiconductors are an important class of materials for 
overcoming the limitations of dilute magnetic semiconductors for spintronics applications. CrSiTe 3 
is a particularly interesting member of this class, since it can likely be exfoliated down to single 
layers, where T c is predicted to increase dramatically. Establishing the nature of the magnetism 
in the bulk is a necessary precursor to understanding the magnetic behavior in thin film samples 
and the possible applications of this material. In this work, we use elastic and inelastic neutron 
scattering to measure the magnetic properties of single crystalline CrSiTe 3 . We find that there is a 
very small single ion anisotropy favoring magnetic ordering along the c-axis and that the measured 
spin waves fit well to a model where the moments are only weakly coupled along that direction. 
Finally, we find that both static and dynamic correlations persist within the ab-plane up to at least 
300 K, strong evidence of this material’s two-dimensional characteristics that are relevant for future 
studies on thin film and monolayer samples. 

PACS numbers: 75.30.Ds, 75.40.-s, 75.50.Pp, 85.75.-d 


I. INTRODUCTION 

The isolation of single layers of carbon atoms, in the 
form of graphene, resulted in the first physical realiza¬ 
tion of a true two-dimensional (2D) crystal leading to 
considerable study and revealing new physic^ 1 . The in¬ 
terest in graphene was not only fundamental, but also 
applied; potential device applications, such as spintronic 
devices, of such two-dimensional materials were recog¬ 
nized soon after the initial discovery 2 . Device applica¬ 
tions of graphene are limited by its intrinsic properties, 
for instance the lack of a band gap, and the search for 
2D crystals soon extended to other materials. Spintron¬ 
ics requires semiconductors that possess ferromagnetic 
properties, particularly those that can be manipulated 
into one- or two-dimensional forms, so much work has fo¬ 
cused on dilute magnetic semiconductors 3 ®, which offer 
a large range of magnetic properties. These compounds 
have been the subject of much debate since their proper¬ 
ties are highly dependent on epitaxial growth conditions, 
dopant distributions and even how the dopants are incor¬ 
porated into the overall band structure--. In an effort to 
overcome these limitations, there has been an emergence 
of research into the much less common spintron ic c andi- 
date: intrinsically ferromagnetic semiconductors 7 ®. The 
compound CrSiTe 3 belongs to this class, having an elec¬ 
tronic band gap of 0.38 meV and a ferromagnetic transi¬ 
tion at T c = 33 K. 

Much of the recent work on CrSiTe 3 has been in 
the form of theoretical calculations focused on the 
band structure and semiconducting properties, particu¬ 
larly calculations of the monolayer properties. General¬ 


ized gradient approximation (GGA) calculations on bulk 
CrSiTe 3 have calculated a bulk band gap of 0.6 meV, 
close to the experimental value, with the gap formed 
by the splitting of the Cr and Te level^. These cal¬ 
culations correctly predict a ferromagnetic ground state 
that induces a slight spin polarization in the Si and 
Te atomiM When the spins were assumed to be Ising- 
like and aligned along the c-axis, local density approxi¬ 
mation (LDA) calculations predicted a nearest-neighbor 
exchange J = -0.58 meV and transition temperature 
T<y = 23 KP®. Experimentally, the structure has been 
determined in earlier work 11 . Hexagonal planes of Cr 3+ 
(S = 3/2) atoms are stacked along the c-axis, with each 
atom octahedrally-coordinated by Te. The Te-Te bond 
lengths are « 3.15 A in the a6-plane and ss 3.48 A out 
of t he pl ane, suggesting a very small octahedral distor- 
tiorPSUl Previous neutron measurements also found a 
magnetic transition at 32.1 K corresponding to the Cr 3+ 
moments ordering ferromagnetically along the c-axiJ^ 3; . 
These measurements found a spin gap of nearly 6 meV, 
which was taken as evidence for Ising behavior. How¬ 
ever, the inelastic measurements suffered from low statis¬ 
tics, which casts doubt on the accuracy of the determined 
magnetic dynamics. Additionally, they do not measure 
in multiple Brillouin Zones, making it possible that they 
have misidentified phonon modes as magnons. This of¬ 
fers an explanation as to why their measurements find 
three spin wave branches, a conclusion that is inconsis¬ 
tent with linear spin wave theory based on the known 
magnetic structure. Additionally, the claim that the 
spins are Ising-like cannot be conclusively derived from 
these measurements for the same reason. Other experi- 
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mental work cited in that work does not contain evidence 
of Isingspins, since the cited magnetostriction measure¬ 
ments 14 are only dependent on the exchange anisotropy, 
while measurem ents o f the single-ion anisotropy show it 
to be quite small 15 Finally, the behavior of the mag¬ 
netic correlations above T c were only discussed qualita¬ 
tively, leaving much room for a more accurate neutron 
scattering study to determine the magnetic correlations. 

The resurgent interest in studying CrSiTe 3 has been 
driven by its applicability to spintronics. This has been 
prompted by the speculation that it may be possible to 
exfoliate the material down to single layeriM which has 
been predicted to have the desirable effect of increasing 
the band gap to « 0.59 e\®and the Curie temperature to 
~ 92 hP 17 ) These features have quickly created the need 
to understand the electronic and magnetic properties of 
CrSiTe 3 . 



II. EXPERIMENTAL DETAILS 


CrSiTe 3 single crystals were grown using a self-flux 
technique, as previously reported 12 . CrSiTe 3 is rhorn- 
bohedral, crystallizing in the space group R- 3, which 
was confirmed with x-ray diffraction measurements on 
the samples used in this study. The crystal structure is 
shown in Fig. [T] Susceptibility measurements showed a 
ferromagnetic transition at T<-< = 33(1) KP^. The large 
c/a ratio and the reasonably accessible magnetic tran¬ 
sition make this a good candidate to study for its ap¬ 
plication as a two-dimensional spintronic material. In 
order to fully characterize the nature of the magnetic 
correlations, including their lower dimensional proper¬ 
ties, neutron scattering measurements were performed 
at the HB-3 and CG-4C (CTAX) triple axis spectrom¬ 
eters of the High-Flux Isotope Reactor, as well as the 
SEQUOIA time-of-flight spectrometer at the Spallation 
Neutron Source of the Oak Ridge National Laboratory. 

Five single crystals of total mass 4.1g and a mosaic of 
2.25° were coaligned in the [H 0 L\ scattering plane for 
use in the CTAX experiment. This array of single crys¬ 
tals was also used for the SEQUOIA experiment, while 
the largest single crystal (mass ~ 1.2g) was used for the 
HB-3 experiment, also aligned in the [H 0 L] scatter¬ 
ing plane. The HB-3 measurements were performed in 
a closed-cycle refrigerator with a base temperature of 
4.0 K using a fixed final energy of 14.7 meV. PG (002) 
monochromator and analyzer crystals were used with PG 
filters, and the collimation was 48’-40’-40’-120’. The SE¬ 
QUOIA measurements were also performed in a closed- 
cycle refrigerator using fixed incident energies of 30 and 
65 meV. The crystals were rotated in the [H 0 L] plane 
in 1° steps over an 83° range. The CTAX measurements 
were performed in a He-4 cryostat with a base tempera¬ 
ture of 1.5 K using fixed final energies of 3 and 5 meV. 
This experiment used a PG002 monochromator, a Be fil¬ 
ter and collimation settings of guide-open-80’-open. 



FIG. 1. (color online) (a) The crystal structure of CrSiTe 3 . 
The Cr 3+ ions form hexagonal arrangements in the a6-plane, 
which are then stacked along the c-axis in an ABC -type 
stacking. This creates two magnetically-inequivalent Cr sites, 
shown in red and blue. The Cr-I site (blue) has a Cr-II atom 
(red) above and a Sb pair (green) below, while it is reversed 
for the Cr-II site. Below T c = 33.2 K, the Cr S = 3/2 
spins align ferromagnetically along the c-axis, due to a single¬ 
ion anisotropy, D,. The nearest-neighbour exchanges out to 
8 A are shown in the figure, which includes 3 in-plane ex¬ 
changes, Jabii Jab 2 and J a b 3 , and one out-of-plane exchange, 
J c . (b) The structure of a single layer is shown, highlighting 
the different exchange interactions in this plane. 


III. ELASTIC NEUTRON SCATTERING 


Fig. §a) shows the temperature dependence of the 
(1 0 1) Bragg peak, which is both a nuclear and a mag¬ 
netic peak, but the magnetic contribution is approxi¬ 
mately 10 times larger at 4 K. The red line is a fit to 
a critical exponent, which yields a transition tempera¬ 
ture of 33.2(1) K and a critical exponent, /3 = 0.151(12). 
This is close to the value expected for a two dimen¬ 
sional transition (housing = 0.125), and well below 
the values expected for a three dimensional transition 
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FIG. 2. (color online) (a) The temperature dependence of 
the measured intensity at the (1 0 1) Bragg peak. This peak 
is both nuclear and magnetic, but the magnetic intensity is 
larger by approximately a factor of 10. A fit to this curve gives 
a transition temperature of T c = 33.2(1) K and a critical 
exponent /3 = 0.151(12). (inset) Scans through the (1 0 1) 
Bragg peak along H at T = 10 K (black squares), 33 K (red 
circles) and 50 K (blue triangles), plotted on a logarithmic 
scale. This shows the sharp increase in the intensity at (1 0 1), 
corresponding to the 3D order, while the diffuse scattering, 
corresponding to the 2D correlations, is most intense at the 
transition, (b) The temperature dependence of the in-plane 
correlation length, determined from the 2-axis scans described 
in the text. The line is a guide to the eye. We note that the in¬ 
plane correlation length is greater than the nearest-neighbor 
Cr-Cr distance at all temperatures measured, suggesting that 
the in-plane magnetic correlations are important well above 
the bulk ordering temperature. 


(/^3 D,Ising — 0.326 and /?3X ),Heisenberg — 0.367). The low 
value of critical exponent is likely a consequence of strong 
two-dimensional correlations in this material. The long- 
range order below 33.2 K is, however, three-dimensional 
in nature, with the spins aligned ferromagnetically along 
the c-axis direction. This was confirmed by measuring 
39 magnetic peaks in the [H 0 L] scattering plane, whose 
intensities below T c were consistent with c-axis ordering. 
This is in agreement with previous work 13 . 

In addition to the magnetic diffraction peaks that 


emerge below Tc, diffuse scattering develops around the 
Bragg peaks, with an intensity that has a maximum at 
T c, shown in the inset to Fig. |2ja). This diffuse scat¬ 
tering arises due to the in-plane correlations that are 
present, particularly above the ordering temperature. In 
order to study the magnetic correlations within the ab- 
plane as a function of temperature, a series of two-axis 
measurements were performed on the HB-3 instrument. 
This was done by measuring along H about the wavevec- 
tor (1 0 0.545), chosen such that the c- axis was parallel to 
the final neutron wavevector, fc/ 18 19i . The analyzer was 
then removed to integrate over various momentum trans¬ 
fers, but the c || kf orientation means that the component 
of the momentum transfer that varies with final energy is 
constrained to be along the L-direction. This provides an 
accurate measure of the energy-integrated spin response 
for two-dimensional correlations within the aft-plane. 

The two-axis measurements are peaked at IT = 1 for 
all temperatures measured and we find that the peak in¬ 
tensity is largest at Tc, which is shown in the inset to 
Fig. [2jb). The width of the peak can be used to cal¬ 
culate the in-plane correlation length, which is plotted 
in Fig. Hfc). As expected, the correlation length peaks 
at T c, but decays as the temperature increases. This 
experiment did not have the required temperature sta¬ 
bility or density of data points near Tp to perform an 
accurate measurement of the critical exponent, v, or to 
observe the expected 2D to 3D crossover, which typically 
occurs within 0.01% to 0.1% of Tc®. However, we do ob¬ 
serve that the correlation length remains larger than the 
nearest-neighbor distance at all temperatures measured, 
up to 250 K. This suggests that short-range correlations 
between the moments in the aft-plane exist well above the 
ordering transition. 


IV. SPIN WAVE MEASUREMENTS 

Fig-i a) and (c) show the spin wave dispersions mea¬ 
sured on the SEQUOIA spectrometer, with an incident 
energy E* = 30 meV. Panel (c) shows the dispersion along 
(0 0 L) with integration ranges of 0.4 reciprocal lattice 
units (r.l.u) along H and K. This shows that the dis¬ 
persion is very weak along the c-axis, and that the mini¬ 
mum in the spin wave dispersion is below 1 meV at the 
zone center, (0 0 -3). Since the dispersion along L is 
relatively weak, it was observed that varying the inte¬ 
gration range along L had little effect on the plot of the 
spin wave spectrum. Therefore, the spectrum shown in 
Fig. |3](a) was constructed by integrating over the full L- 
range, —12 < L < 12, while using an integration range 
of 0.2 r.l.u. along K. This shows two spin wave bands, 
one peaking at ss 8 meV and the higher band peaking at 
ss 15 meV, which meet at the zone boundaries. The data 
along {H 0 0) using a smaller integration range for L 
showed identical dispersions, supporting the conclusion 
that the spin waves are close to two-dimensional. The 
spin waves were modeled with a Hamiltonian given by: 
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FIG. 3. (color online) The panels on the left show the spin waves of CrSiTe 3 measured on the SEQUOIA spectrometer at 
T = 10 K, while the panels on the right show the calculated spin wave dispersions, along the (H 0 0) direction (panels (a) and 
(b)) and the (0 0 L ) direction (panels (c) and (d)). The calculated patterns use the model and exchange parameters described 
in the text, producing very good agreement with the experimental data. 


H — Jabl ^ ' Si * jS )-|-1 T J a b2 ^ ' Si * Si+2 

i i 

+ Jc^Si. S l+ 3 + J ab3 Y Si ■ s i+ 4 -D Z Y (S -) 2 

( 1 ) 

where the J’s are exchange constants between neighbor¬ 
ing spins, Si, D z represents a single-ion anisotropy and 
the spins are assumed to be localized Cr 3+ , S = 3/2 mo¬ 
ments, as shown in Fig. |Tj The very small value of the 
spin gap and the weak c-axis dispersion are both able to 
be accurately determined from the SEQUOIA data. The 
magnetic anisotropy necessitates a non-zero spin gap, but 
the SEQUOIA measurements have shown that it is less 


than 1 meV. This is consistent with magnetization data 
taken to characterize the sample, which observes very lit¬ 
tle magnetic hysteresis and a field of 1.5 T applied in the 
afr-plane (the hard magnetization direction) is required to 
saturate the magnetic moment. Additionally, DFT calcu¬ 
lations based on purely van der Waals interactions have 
found that there is an entropy difference of « 20 iteV, 
favoring c-a xis o rdering, as compared to ordering within 
the a6-plan^- 21 i This is very close to the value obtained 
for the single-ion anisotropy, as discussed below. The 
fits described below produce the simulations shown in 
Fig. |3|b) and (d). These show very good agreement with 
the experimental data, but we do observe a discrepancy 
at the non-zero integer values of H. At these points in the 
experimental data, we observe broadening of the optical 
mode and intensity shifts to the acoustic mode. These 
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FIG. 4. (color online) (a) The constant-I? scans measured on the Cold Triple-Axis (CTAX) Spectrometer at 1.5 K along 
(0 0 L), which are offset along the y-axis for clarity. These scans focus on the region near the zone center (0 0 3), where 
a spin gap is expected in the spectrum. By assuming that the energy of the spin waves, E oc |q| 2 in the long wavelength 
limit, we can extract an approximate measure of the gap, despite it being within the resolution of the instrument, (b) The 
q-dependence of the spin waves near the zone center is plotted, showing g 2 -dependence. A fit to these values gives a value for 
the gap A = 0.075(24) meV. (c) A constant-Q measurement at Q = ( 1.7 0 0) near the zone boundary, (0 0 1.5), allows for a 
determination of the value of the out-of-plane coupling, J c = -0.730(96) meV. 


features are not seen in the simulation, though the origin 
of this difference is unclear. 

In order to measure the magnitude of the spin gap and 
the dispersion out of the plane, measurements were per¬ 
formed on the cold neutron triple axis (CTAX) instru¬ 
ment at Oak Ridge National Laboratory. The energy 
resolution of a constant-Q scan was insufficient to sep¬ 
arate the gapped mode from the elastic scattering. To 
extract the spin gap value, a series of constant-i? scans 
was performed, shown in Fig. |4|a). These scans focus 
on the region near the zone center (0 0 3), where a spin 
gap is expected in the spectrum. By assuming that the 
energy of the spin waves, E oc |< 7 |“ in the long wavelength 
limit, we can extract an approximate measure of the gap, 
despite it being within the resolution of the instrument. 
The data was fit in multiple zones, as shown in Fig[4](b), 
showing g 2 -dependence. This gave a value of the spin 
gap, A = 0.075(24) meV, allowing a measurement of the 
single-ion anisotropy, D z = A/2 S = 0.0252(80) meV. 
The single-ion anisotropy arises from the small c-axis dis¬ 
tortion in the Te octahedra surrounding the Cr 3+ ions. 
This nearly perfect octahedral environment makes spin- 
orbit effects very weak and the Cr-Cr bond is a point 
of local centrosymmetry, ruling out Dyaloshinskii-Moriya 
interactions as an origin for the energy gap. This small 
gap is consistent with theoretical predictions that the 
energy difference between c-axis and a6-plane ordering is 
less than 0.1 meV. Constant-Q measurements near the 
zone boundary, (0 0 1.5), shown in Fig. [4j)c), allow us 
to determine the value of the out-of-plane coupling, J c , 
since the energy of the lower spin wave branch along L 
depends only on D z and J c . Using the value of D z de¬ 
termined above, we obtain J c = -0.730(96) meV. 

These values were then fixed for the purposes of fitting 
the SEQUOIA data to obtain the in-plane exchange con¬ 
stants. The fitting was performed using the fitting rou¬ 
tines built into the Horace software package^ which used 


the dispersions calculated by SpinW^as its model, com¬ 
paring the intensities of the data and calculation over the 
entire range in the Q-E slice. This allowed the in-plane 
exchange constants to be fit simultaneously, utilizing the 
two-dimensional dataset shown in Fig. [3^a) . The values 
of these three exchange constants obtained from fitting 
the data shown in Fig. [3](a) to Eq. [lj as well as those 
determined from the CTAX measurements, are given in 
Table H 


Exchange 

Description 

Distance 

Value (meV) 

Jabl 

1 st NN in-plane 

3.907 A 

-1.27(23) 

Jabl 

2 nd NN in-plane 

6.768 A 

-0.10(50) 

Jab3 

3 rd NN in-plane 

7.814 A 

-0.285(73) 

Jc 

1 st NN out-of-plane 

6.852 A 

-0.730(96) 

D z 

Single Ion Anisotropy 


0.0252(80) 


TABLE I. The exchange constants obtained from fitting the 
inelastic measurements using Eq. [l] To properly describe the 
spin waves, it was necessary to use couplings up to 8 A, which 
requires 3 in-plane interactions and one out-of-plane. Consis¬ 
tent with the quasi-2D nature of the material, the value of 
J c gives a very small dispersion along the c-axis and, as ex¬ 
pected from the very small octahedral distortion, the value of 
the single-ion anisotropy, D z , is very close to zero. 

As expected for a quasi-two-dimensional system, the 
nearest neighbor in-plane coupling, J a bi, is the domi¬ 
nant interaction and is ferromagnetic, as are the other 
in-plane exchange constants. While J a b\ is only slightly 
larger than J c , there are 3 nearest neighbors in the ab- 
plane compared to only one along the c-axis. This makes 
the in-plane exchange more than 5 times larger than 
the out-of-plane coupling. The ferromagnetic origin of 
the in-plane terms is likely due to superexchange medi¬ 
ated by the Te ions, shown in Fig. ITVb) . The Cr 3+ ion 
direct exchange is antiferromagnetic, though the large 
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FIG. 5. (color online) The temperature dependence of spin dynamics, (a) The inelastic neutron spectrum measured at 
T = 40 K on the SEQUOIA spectrometer. We see that broadened spin wave features remain, even above the transition, while 
there are no such excitations along the L direction. This is a clear signature of the two-dimensionality of this material. The 
feature at 20 meV is extrinsic to the sample, verified by performing an empty can measurement, (b) The in-plane spin dynamics 
along (H 0 3) at various temperatures, measured on CTAX and offset for clarity. Sharp spin waves are observed below the 
transition, but dynamic correlations exist at all temperatures measured. The high temperature features are reminiscent of the 
spin waves seen below Tc, but significantly broadened. The lines are guides to the eye. (c) In contrast, there are no dynamic 
correlations along the L direction above Tc (also shown offset for clarity). This suggests that three-dimensional correlations 
only exist below the Curie temperature, while the in-plane dynamic correlations persist up to at least 300 K. 


nearest-neighbor Cr-Cr distance of 3.987 A should make 
this a weak effecP^HD Finally, the Cr octahedra are 
edge-sharing and the Cr-Te-Cr angle is 88.8°, very close 
to perfectly orthogonal. This arrangement of the ions 
would suggest that superexchange of the type described 
by the Goodenough-Kanamori ruletP^HSU would be ferro¬ 
magnetic, consistent with the observations of ferromag¬ 
netism in other layered Cr 3+ compounds with similar su¬ 
perexchange pathwayd^H. First principles DFT calcula¬ 
tions on monolayer CrSiTe 3 suggest that the 2 nd and 3 rd 
nearest-neighbor in-plane interactions are Cr-Te-Te-Cr 
double-superexchange interactional, leading to the small 
values we observe for J a h 2 and J a b3- Furthermore, the ob¬ 
servation that | Jab 2 1 < | Jab 3 | was also predicted by the 
DFT calculations. There are two double-superexchange 
pathways that contribute to J a b2 , one of which is ferro¬ 
magnetic and the other antiferromagnetic, making it a 
very weak exchange 17 . 

The neutron scattering measurements are also consis¬ 
tent with the theoretical predictions of an enhanced tran¬ 
sition temperature in CrG eTes and CrSnTe 3 , as well as 
monolayer CrSiTeJ^ 1 '^. Despite the spins being only 
weakly Ising-like, the easy-axis spin anisotropy created 
by the imperfect Te octahedra still allows for ordering 
in two-dimensions, such as when the materials are re¬ 
duced to monolayers. The loss of the out-of plane cou¬ 
pling is compensated by an increased in-plane ferromag¬ 
netic exchange, as predicted by several theoretical studies 
that compare bulk CrSiTe 3 to its monolayer version and 
the chemically-substituted CrGeTe 3 and CrSnTej- 8 | 17 l 21 l. 
The increased ferromagnetic exchange within the ab- 
plane arises due to the increase of the nearest-neighbor 
Cr-Cr distance in all of these cases. This serves to re¬ 
duce the Cr-Cr direct exchange, which is antiferromag¬ 
netic, while pushing the Cr-Te-Cr bond angle closer to 


90°. This increases the effect of the ferromagnetic su¬ 
perexchange, resulting in the increased values of T C m 
Inelastic neutron scattering measurements performed 
above the ferromagnetic transition are shown in Fig. [5] 
Measurements at 40 K along the (H 0 0) direction are 
presented in Fig. [5](a). This data indicates that the spin 
waves are still present within the planes, but have been 
significantly broadened. The minima in the spin waves 
seen at (1 0 0) and (2 0 0) in Fig. [3](a) have disap¬ 
peared, as they arose from the weak dispersion along the 
T-direction. This indicates that there is no dispersion 
along L above Tc, and these broadened spin waves are 
only present within the afr-planes at these temperatures. 
To track the temperature evolution of the spin waves, 
constant-T measurements were performed along {H 0 0) 
at E = 1.75 meV and (0 0 L) at E = 1.5 meV as a 
function of temperature. These are shown in Fig. [5jb) 
and (c), respectively. We see that the spin waves along 
L disappear abruptly above Tc, but that the broadened 
spin waves exist at the same Q along (H 0 0) well above 
the ordering transition, up to at least 300 K. 


V. CONCLUSIONS 

We have measured CrSiTe 3 using elastic and inelastic 
neutron scattering. These measurements observe bulk 
ferromagnetism below Tc = 33.2(1) K and the critical 
exponent of the neutron scattering intensity was found 
to be /3 = 0.151(2), close to the value expected for a 
two-dimensional system. The magnetic Bragg peaks also 
exhibited a diffuse component, the intensity of which 
peaks at Tc, indicating two-dimensional ferromagnetic 
correlations are present in the a&-plane above the ferro¬ 
magnetic transition. To characterize these correlations, 
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two-axis measurements were performed, which provided 
an effective quantitative measure of the in-plane correla¬ 
tion length. As expected, the correlation length diverged 
at Tc but was still larger than the nearest-neighbor Cr- 
Cr distance at all temperatures measured, up to 250 K. 
These measurements suggest that while CrSiTe 3 only or¬ 
ders in three dimensions below 33 K, there are strong 
two-dimensional static correlations that persist up to at 
least room temperature. 

In contrast to the previous assumption of Ising spins, 
the spin wave measurements suggest that the spins are 
very Heisenberg-like, with a very small spin gap of 
0.075(24) meV due to a small single-ion anisotropy. X- 
ray diffraction measurementf 12 and DFT calculation^ 21 ! 
find evidence of van der Waals interactions creating a 
very small octahedral distortion along the c-axis, creat¬ 
ing the anisotropy. We find that the dominant inter¬ 
action, i = -1.27(23) meV, is likely due to the su¬ 
perexchange mediated by the Te ions. The 2 nd and 3 rd 
nearest-neighbor in-plane interactions are much weaker, 
likely a result of double-superexchange interactions. 

Above Tc, no spin waves are present along the L di¬ 
rection, but they persist in a broadened form along the H 
direction. Measurements at 40 K along (H 0 0) showed 
broadened spin waves consistent with dynamics confined 
to the afo-planes, a clear indication of two-dimensional 
behavior, while temperature-dependent measurements 
show an inelastic signal that is visible up to 300 K. The 
latter finding suggests that there are dynamic in-plane 
correlations that persist up to at least room temperature 
as well. Considering these results in conjunction with 
the neutron diffraction measurements, we find that both 


static and dynamic in-plane magnetic correlations exist 
up to at least 300 K, ten times Tc- This illustrates the 
importance of two-dimensional correlations in these ma¬ 
terials and may suggest that similar physics is relevant 
in the isostructural compounds CrGeTe 3 and CrSnTe 3 , 
as well as for monolayer CrSiTe 3 , where it has been pre¬ 
dicted that increased separation or complete decoupling 
of the hexagonal Cr layers puts an increased emphasis on 
the two-dimensional correlations. Taken together with 
the dramatic increase in the magnetic transition temper¬ 
ature driven by a separation of the intralayer Cr atoms 
in these materials, this suggests that spintronic device^ 
that make use of these semiconductors on a substrate 
that enhances the magnetic character through strain or 
other effects may be capable of supporting intrinsic fer¬ 
romagnetism at temperatures approaching, or even sur¬ 
passing, room temperature. 
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